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Abstract: For many practical purposes, it is convenient to formulate unbroken 
non-abelian gauge theories like QCD in a color-flow basis. We present a new deriva- 
tion of SU(iV) interactions in the color- flow basis by extending the gauge group to 
U(iV) x U(l)' in such a way that the two U(l) factors cancel each other. We use the 
quantum action principles to show the equivalence to the usual basis to all orders in 
perturbation theory. We extend the known Feynman rules to exotic color represen- 
tations (e. g. sextets) and discuss practical applications as they occur in automatic 
computation programs. 
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1 Introduction 

The analysis of particle physics experiments at colliders depends on reliable theo- 
retical predictions for cross sections of scattering processes. In the LHC era, hard 
processes with more external partons than ever have become accessible, and their 
analysis is essential for unveiling the physics at the Terascale. This results in a two- 
sided challenge: on one hand, complex Standard Model processes must be computed 
with unprecedented precision and on the other hand, the parameter space for many 
models for physics beyond the Standard Model must be scanned with sufficient ac- 
curacy. This situation has prompted the development of computer programs that 
can automatically compute differential cross sections and sample them efficiently on 
phase space starting from the specification of a Lagrangian or, equivalently, a set of 
Feynman rules. 

In QCD, quarks and antiquaries come in three colors and gluons in eight. Thus, 
there are many amplitudes with different color quantum numbers that must be com- 
puted and summed in quadrature, eventually. In Feynman diagram based calcula- 
tions, the contribution of each diagram factorizes into the dependence on color and 
the dependence on flavor, polarization and momenta. This allows to compute the 
color factors once and for all. However, for processes with many external particles, 
the number of Feynman diagrams grows factorially, and more efficient, e.g. recur- 
sive, algorithms must be used that can take advantage of the cancellations among 
diagrams. Unfortunately, the color dependence does not factorize for the whole 
scattering amplitude and new expansions, like color ordered amplitudes must be 
employed to allow a separate computation of color factors. 

It is therefore worthwhile to investigate efficient representations of the color de- 
pendence that work for complete amplitudes, i. e. sums of Feynman diagrams. It 
turns out that expressing everything in reducible tensor products of the fundamen- 
tal representation and its conjugate instead of higher irreducible representations is 
beneficial both for computation and for interfacing to parton showers, fragmentation 
and hadronization. If there are no interaction vertices with exotic color structures, 
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the resulting amplitudes can always be decomposed into weighted sums of products 
of Kronecker-5s, so called color flows 1 . 

The adjoint representation of U(iV) is isomorphic to the tensor product N x N 
of the fundamental representation and its conjugate. Thus, the color coordinate a of 
a U(A r )-gauge boson in the adjoint representation can be represented equivalently by 
a pair of using the decomposition A/ = A a [T a ] i J . This equivalence is employed 
in 't Hooft's double line notation [1], in which each U(iV)-gauge boson line is replaced 
by a pair of matter and anti-matter lines for the purpose of computing color factors: 



It was soon realized that the diagrammatic rules could be extended to accomodate the 
tracelessness of the generators of SU(iV) [2]. Indeed, an efficient diagrammatical al- 
gorithm for computing arbitrary color factors in any gauge group was constructed [3] 
using a generalization of the double line notation to higher representations 2 . 

However, the color-flow representation is more than a useful calculational de- 
vice: Feynman diagrams in the double line notation bear a striking resemblance to 
dual diagrams [6, 7] and indeed the iV — > oo-limit is related to a theory of open 
strings [1]. It is therefore no accident that the successful semiclassical models for the 
fragmentation of quarks and gluons that are based on string dynamics [8] require 
a specification of the parton cross sections in the color-flow basis [9]. The earliest 
examples of Monte-Carlo event generators that combined the partonic cross sections 
for the hard interactions with QCD-inspired fragmentation models [10] provided a 
template for the interfacing of hard interactions with fragmentation. Similarly, a 
color-flow representation is used in leading order in 1/N for the implementation of 
color coherence in parton showers [11]. These interfaces became the standard HEPEVT 
common block [12], which includes color flow information. Originally developed for 
interfacing Monte-Carlo event generators for hard scattering matrix elements with li- 

x If there are couplings like a totally antisymmetric vertex Y2i j k=i ei 3k$i < f > j<f > k m SU(3), the 
corresponding invariant tensors can also appear in the final result. 

independently of Feynman diagrams and quantum field theory, Penrose had earlier introduced 
a diagrammatical notation for general tensor calculus, including group representations [4]. This 
approach has evolved into the backbone of an unconventional textbook on the exceptional Lie 
algebras [5]. 
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braries or external programs for parton showers, fragmentation and hadronization at 
LEP1, it was widely adopted and has grown into the Les Houches Accords (LHA) [14]. 

Therefore the color-flow basis is particularly suited for the computation of par- 
tonic scattering amplitudes that will be interfaced with the other components of a 
Monte-Carlo event generator. Indeed, most generators for hard multi-parton scat- 
tering can provide the corresponding information. One can, however, go one step 
further and not only represent the final result in a color-flow basis, but use Feynman 
rules in a color-flow basis already for the computation. A lucent derivation of these 
Feynman rules was given in [15]. Since the derivation in [15] is based on the QCD 
Lagragian, it is not necessary to use Feynman rules at all, but one can implement 
modern recursive algorithms [16-18] directly in the color-flow basis [19-21]. 

However, the derivation [15] of the color-flow representation is incomplete in two 
directions: firstly, there is no consideration of interactions with more exotic color 
structures, in particular beyond QCD with fermionic matter in the fundamental 
representation, and secondly the discussion is deliberately confined to tree level am- 
plitudes. These two limitations are in fact related, and overcoming them is not only 
of theoretical interest: the most important light Higgs production channel at LHC 
involves the dimension-5 operator H Tr (F^F^), which arises from a loop and corre- 
sponds to an octet-octet-singlet coupling that cannot be described straightforwardly 
in the framework provided by [15]. Indeed, our generalization of the color-flow rep- 
resentation was prompted by the implementation of the effective Higgs-gluon-gluon 
vertex in WHIZARD's [19] optimized matrix element generator O'Mega [18]. While 
it was simple to "fudge" the Feynman rules for a single insertion, it turned out that 
these prescriptions would give incorrect results for multiple insertions. 

This paper is organized as follows: we start with a motivation and descrip- 
tion of our implementation of color-flow QCD in section 2 and discuss quantum 
field-theoretical aspects in more detail in section 3. In section 4, we show that our 
description is valid to all orders in perturbation theory. In the following two sections 
we discuss some applications in detail: tree level amplitudes in section 5 and the 
already mentioned effective interactions from loops in section 6. For reference, we 
repeat the color-flow Feynman rules for QCD with fermionic matter in the funda- 
mental representation in appendix A. As an example of exotic color representations, 
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we extend those rules to color-sextet particles in appendix B. 
2 The Model 

In this section, we construct a variant of QCD which contains two extra gauge degrees 
of freedom. The model is set up not for its physics content, but for its usefulness 
in practical perturbative calculations. In fact, as we will show in later sections, in 
physical quantities the extra degrees of freedom cancel, so all predictions are identical 
to ordinary QCD. 

2.1 QCD Preliminaries 

We consider QCD with a single fermionic matter species in the iV-dimensional fun- 
damental representation of SU(iV). While we are really interested in the N = 3 case, 
it is sometimes useful to explicitly keep the dependence on N, as we will do through 
most of this paper. Furthermore, we will not assume that the matter representa- 
tion is vector-like as in QCD, although we do follow the QCD notation with Dirac 
fermions in the Lagrangian. 

Let us first recall basic facts of perturbative QCD. The perturbation series is 
derived from a Lagrangian which splits into two parts, 

£ = Anv + £gf, (2.1) 

a gauge-invariant part £ inv and a gauge-fixing part £ g f. The gauge-invariant part is 
given by 3 

A™ = -^Tr G^& v + 4, + 4) if>. (2.2) 

The gauge-fixing part depends on the chosen gauge-fixing procedure. In a manifestly 
covariant formulation with linear gauge-fixing, it takes the form 

C g{ = \TrB(d ■ A) + 4rTr B 2 - ^Tr td^c - i[i M , £]). (2.3) 

Since the matter fields are in the fundamental (defining) iV- dimensional representa- 
tion of SU(iV), we choose to represent the Lie-algebra valued fields G, A, B, c, c by 
traceless N x N matrices. 4 The field-strength tensor G can be expressed in terms of 

3 We set the masses of matter fields to zero, for brevity. 

4 For later convenience, we mark traceless matrix fields by a hat, A v . 
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the gauge potential A as 

G>, = d„A u - d v A» - i[A„, A v \. (2.4) 

Furthermore, the gauge-fixing term involves a Fadeev-Popov [22] ghost field c, an 
antighost field c and a Nakanishi-Lautrup (NL) [23] auxiliary field B. The latter 
can be integrated out in order to obtain the more familiar form of the gauge fixing 
Lagrangian 



£ gf = -— Tr (d ■ A) 2 - -Tr S9"(fyc - i[A M , c]). (2.5) 



^Tr(9-i) 2 -4 

We may assume a manifestly gauge-invariant renormalization procedure such as 
MS, so that the Lagrangian retains the form (2.2,2.3) in each order of the perturbative 
expansion. 5 The fields in the Lagrangian are understood to be renormalized, order 
by order, such that the model contains only two (renormalized) real parameters, the 
gauge-coupling g and the gauge-fixing parameter £. The latter drops out of physical 
quantities. 

The gauge group SU(iV) is a subgroup of the general linear group GL(iV), there- 
fore the fields (e.g., A), in the N x N matrix representation, obey algebraic con- 
straints: 

A ] = A and Tri = 0. (2.6) 

The constraints are automatically satisfied if we introduce the usual basis for the Lie 
algebra representation, T a (a = 1, . . . N 2 — 1), and write A = ^2 a A a T a . The basis 
elements obey the hermiticity and trace conditions: 

Tr T a = 0, Tr T a T b = \5 a \ [T a , T b ] = if abc T c , (T a ) t = T a (2.7) 

with the real structure constants f abc . Each Feynman graph can be factorized into a 
kinematical and a color amplitude. The color amplitude consists of a string of T and 
/ tensors, contracted over all internal indices, and representable as a combination 
of tensors with open external indices a,b, . . . ,i,j, ■ ■ ■■ For squared amplitudes, all 
indices are summed over, hence the particular representation becomes irrelevant. 
Various schemes exist to compute color amplitudes exactly and efficiently. 



5 A generic renormalization procedure could make non-invariant terms in the renormalized La- 
grangian necessary that ensure gauge invariance of the effective action. 
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2.2 The Color-Flow Representation 

In the following, we will eliminate the basis T a , and instead work with the individual 
matrix elements of the gauge fields in the chosen Lie algebra representation [5]. This 
approach is known as a color-flow representation. For concreteness, we consider the 
gauge potential A. Without algebraic constraints, a N x N matrix field contains N 2 
complex degrees of freedom. 

We have to implement the algebraic contraints (2.6). The hermiticity condition 

is 

&3 = (2.8) 

In effect, each unordered index combination (ij) carries only one, instead of two, 
complex degree of freedom. We can understand this as if each ordered index com- 
bination carried one real degree of freedom. For the propagator in a Feynman 
graph, a value of an index can be represented by a colored line. There are N different 
colors which correspond to the N values each index may take. An arrow of the line 
indicates whether the index is the lower or upper one. Hence, A propagators carry 
two color lines with opposing directions. 

Similarly, matter propagators are represented by a single directed color line. 
External states are represented by a terminating line (matter) or two opposing ter- 
minating lines (gauge). As long as we consider only the interactions following from 
the Lagrangian as quoted above, the rule of matrix multiplication trivially ensures 
that color is conserved at each vertex, i.e., lines end only at external states. 

While this graphical approach is straightforward, we also have to enforce the 
trace condition from (2.6), explicitly 

£i\ = 0, (2.9) 

which reduces the number of degrees of freedom by one. The implementation of this 
constraint amounts to additional modified color-flow patterns in Feynman graphs [15]. 
While the construction of these patterns is easy to understand, they are, by definition, 
a non-local modification of the naive color-flow expansion. This is a complication in 
automatic, in particular non-diagrammatic, algorithms, which we want to avoid. 

Instead, we follow a field-theoretical approach and set up a modified QCD theory. 
This theory has the algebraic constraints inherent in its field content, such that they 
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are automatically satisfied by the generated color-flow amplitudes. In essence, it 
generates (i) terms that complete the A field to a U(iV) matrix so that it has no 
trace condition and naive color-flow Feynman rules apply, and separately (ii) terms 
that subtract the spurious contributions without complicating the color-flow rules. 
To make the diagrammatic expansion unambiguous, the two fields have to be formally 
independent of each other. 

2.3 Singlet and Phantom Gluons 

We denote the field which complements the gluon field matrix as the singlet gluon A . 
This field is associated with a "zeroth" generator T° which is not traceless but satisfies 



TrT°=Jf, TrT a T° = ±5 a0 , [T a ,T°] = 0, (2.10) 
and let a run from to N 2 — 1. Explicitly, 



(T°V. = (2.11) 




(We choose upper indices for the fundamental and lower indices for the anti-fundamental 
color representation.) The matrix field resulting from combining A with A Q 



(2.12) 



has N independent components, the color- flow gluons A 1 y Analogously, for the NL 
and ghost fields, we replace the traceless fields by unconstrained fields B,c,c and 
singlet fields Bs, cs, cs, respectively. The field strength G is modified accordingly. 

For the subtraction terms, we introduce a phantom gluon field A. The phantom 
gluon is an independent U(l) gauge boson with the wrong sign in the propagator. 
Like the other gluons, it couples to the matter fields via the QCD coupling g. It does 
not couple to gluon fields. We also introduce the corresponding NL field B. (We 
may also introduce phantom ghost fields, but these can be omitted with linear gauge 
fixing, since the U(l) gauge group is abelian.) 

2.4 Color-Flow QCD 

We can now set up the Lagrangian for the color-flow version of QCD, that includes 
the singlet gluon (as a component of the gluon matrix A* -) and the phantom gluon 
A, together with their associated field strength, NL and ghost fields. 
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The Lagrangian 6 is 

C = -^Tr G^W + ^fi^GV + i> + A - A) 1> 

2 2 AT £ £N ~ 

+ -Tr B(d ■ A) - —B(d ■ A) + ^Tr B 2 - ^—B 2 + £ ghost (2.13) 

g 2 g 2 g 2 g 2 

where 

^ghost = -^Tr cd^c - i[A„, c]) (2.14) 

and the decoupling ghosts for the phantom gluon have been omitted. 

A priori, this theory is different from QCD. We will denote it as "color-flow 
QCD" and we will show its equivalence to QCD in section 3. It contains N 2 instead 
of N 2 — 1 independent degrees of freedom in the A field, plus the additional U(l)' 
gauge field A with wrong-sign propagator. Thus, it can be regarded as a U(iV) x U(l)' 
or, equivalently, SU(iV) x U(l) x U(l)' gauge theory. In a path-integral formulation, 
there are independent integrations over all N 2 + 1 gauge components. The singlet 
field A is present, but hidden in the gauge- field matrix. 

In the Lagrangian (2.13), we have encoded the algebraic properties of the color- 
flow decomposition as a (perturbative) quantum field theory on its own. Therefore, 
we can expect that amplitudes calculated from this field theory manifestly exhibit 
the advantages of this approach. 

In calculations, we therefore retain the decomposition of the matrix-valued fields 
in terms of their matrix elements, A 1 -. As explained above, these non-hermitean vec- 
tor bosons are represented by double lines, one for the color and one for the anticolor 
flow, flowing in opposite directions. The vertex color factors are Kronecker deltas, 
which diagrammatically correspond to color lines continuing through the vertex. The 
phantom gluon, on the other hand, is analogous to a photon. In particular, it carries 
no color flow. Color lines always terminate at external states. Hence, for any tree 
diagram, the color factor is simply N p , where p is the number of distinct color lines 
in the diagram (after squaring). Since this is determined by the colors of the external 

6 For the purpose of discussing renormalization, we incorporate the gauge coupling in the nor- 
malization of the gauge-field and gauge-fixing terms. To derive Fcynman rules later, we will apply 
a trivial field renormalization to transfer g/V2 to the vertices and obtain a canonically normalized 
U(A r ) gauge field from (2.13). However, we will retain the factor N in the phantom kinetic term 
and it will appear in the phantom propagator and external states. 
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states, the whole color algebra becomes trivial. The price for this is a proliferation 
of extra Feynman diagrams (which have identical kinematics, however, and thus do 
not blow up the calculation.) 

We list the Feynman rules that follow from (2.13) in Appendix A. 

2.5 Example 

Let us look at the exchange of a gluon between fermion lines. In ordinary QCD, this 
is represented by the diagram 




In color-flow QCD, we get two diagrams, one for gluon and one for phantom exchange 




with the correct sum 



' a 

Usually, this color- flow decomposition is obtained by applying the relation (2.15) 
backwards. The Feynman rules derived from the color-flow Lagrangian generate it 
directly. 

Also at the squared- amplitude level, the color-flow representation leads to a reor- 
ganization of the color-flow calculation, which considerably simplifies the algorithm. 
There is a further benefit: common parton-shower and hadronization algorithms ex- 
pect the color connections to be expressed in a color-flow basis. We obtain both the 
exact amplitude (at the given order) and, by omitting the interference terms and 
all terms containing external phantoms, the projection on the possible color-flow 
patterns in a single step. 
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In summary, we work with a trivial version of color algebra since the complete 
color flow information is represented on a graphical basis. Subtractions in the color 
factors are accounted for by extra diagrams which involve the phantom gluon. No 
approximation is involved. 

3 Field-theoretical Considerations 

In the following, we demonstrate that color-flow QCD as described above is equivalent 
to the original SU(iV) theory, i. e., A and A cancel each other in physical quantities. 
This equivalence may appear trivial, at least in tree-level amplitudes. However, it 
bears some subtleties, in particular if it is to be extended to all orders. To promote 
the approach from a technical device to a consistent field theory, we have to prove 
that the equivalence can be maintained, by suitable renormalization conditions, to 
all orders in perturbation theory. 

This situation is similar to the covariant quantization of QED, where the scalar 
and longitudinal photon components, which are both unphysical, cancel in physical 
quantities. In a nonabelian theory, this requires the additional introduction of ghost 
fields [22, 24-26]. In the Fock space, the unphysical fields lead to indefinite- norm 
states, but it is possible to build a Hilbert space out of equivalence classes of physical 
states in which the 5" matrix is unitary. As we will show below, the same argument 
applies in the present case: the phantom gluon generates negative-norm states, but 
it is always produced coherently with the singlet gluon (which is now hidden in the 
N x N gluon matrix), and therefore a unitary S matrix is well defined. 

3.1 Algebraic Identities at Lowest Order 

A consistent formulation of a quantum field theory is based on the symmetries that 
the effective action, including the integrated Lagrangian as its tree approximation, 
has to satisfy order by order in a perturbative expansion. These are representable as 
algebraic identities [26]. 

In the present subsection, we look only at the lowest order, so we temporarily 
ignore the complications that necessitate the replacement of gauge symmetry by 
BRST symmetry, ghost fields, etc. Given the fields A l j, A, B^, B, ip, ip, we observe 
that the lowest-order action T° = J C from (2.13), with the ghost term omitted, 
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satisfies the Ward identity 7 of U(iV) 
Tr 

where the gauge transformations of the fields are 



— Tr [Bd ■ S A] 



(3.1) 



SA = du-i[A, u], SA = 0, Sip = iuip , Sip = -iipu , SB = . (3.2) 

This Ward identity decomposes into separate identities for the SU(iV) and U(l) 
factors. Even though we are interested only in the U(l) part, we want to keep our 
notation free of the clutter of indices. Therefore, we retain the matrix notation and 
use traces to project on the U(l). Indeed, we can evaluate the Ward identity (3.1) 
for ui = ojol to find 

2 



-d 2 TrB. 



This is accompanied by the gauge-fixing condition for the singlet 

Tr — = — (Tr d ■ A + £Tr B) . 

SB g z 



(3.3) 



(3.4) 



Furthermore we define, for any matrix- valued field M G {A, G, . . .}, the linear com- 
binations 



M± = -= ( — TrM± M 
V2\N 



and the dual functional operators 

ST 



ST , ST 



f Tr - 



with 



SM± y/2 V 5M SM 



smx) -f\*-v), SM±{X) 



0. 



SM ± (y) v y " SM T (y) 
Using (3.5), we can rewrite terms appearing in the Lagrangian as 

-Tr M 2 - —M- = -Tr M 2 + NM + M_ 
2 2 2 + 

T°M° - M = V2M_. 



(3.5a) 

(3.5b) 
(3.6) 

(3.7a) 
(3.7b) 



7 A11 functional derivatives are defined to act from the left which introduces a sign for the variation 
of fermions 

OWa Olpa dip Olp 
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Computing the derivative of the lowest-order action with respect to A and A 
and combining the results, we obtain the phantom equation 

* (d„G^ - Nd v BS) . (3.8) 



SA +V v " 

It expresses the cancellation between singlet and phantom terms. Since A + is an 
abelian gauge field, applying 5/b~A v + to the ghost term (2.14) yields zero, so the 
phantom equation holds for the complete lowest-order Lagrangian in the BRST for- 
malism [24]. 

3.2 Fock-Space Cancellation 

To clarify the physical implications of the phantom equation (3.8), let us use (3.7) 
to express the singlet and phantom fields (^Tr A, A) by (A + , AS). The remainder of 
the matrix A is traceless, this is the original SU(iV) gluon with all of its interactions. 
We must show that the extra fields A + and A_ do not contribute to observable 
quantities. 

We eliminate the B fields via their equations of motion. The terms in the re- 
sulting Lagrangian that depend on A + or A_ take the form 

£ = ~G+^G.^ - ^(d ■ A+){d ■ A.) + V2$ (A_) ^ . (3.9) 

The propagator interchanges A + and A_, but only A_ couples to matter. Couplings 
to ordinary gluons are absent in (3.9), because G± = d^A^. — d u A±. As a result, 
whenever a singlet /phantom gluon is created by a matter current, it cannot be an- 
nihilated, and vice versa. In short, A + and A_ do not introduce any observable 
interactions. 

In fact, we may freely add any term to the effective action that depends on A^ at 
least linearly, but has no dependence on A + . Any such interaction is unobservable. 

This pattern will continue to hold to all orders, unless a loop diagram induces an 
interaction of some current with A + . However, the phantom equation (3.8) prohibits 
this: it tells that the interactions of A + are completely accounted for by the leading- 
order Lagrangian, up to a renormalization of the coupling constant. Therefore, if we 
can maintain the phantom equation on the effective action to all orders, observables 
computed from the color-flow QCD theory are identical to those from ordinary QCD. 



13 



In actual applications, we do not express the Lagrangian in terms of A + and 
A_. Instead, we work with the U(iV) field A and the U(l)' phantom A, as we did 
in the previous sections. Hence, the Fock-space cancellation is not entirely trivial 
since it involves cancellations among graphs that belong to different gauge groups, 
at least superficially. We have to enforce the phantom equation explicitly in order 
to guarantee that no A + interaction arises in the effective action. 

3.3 Symmetries 

If the phantom equation (3.8) is combined with the Ward identity (3.3) and the 
gauge-fixing condition (3.4), we obtain an analogous Ward identity and gauge-fixing 
condition for the phantom gluon. The complete symmetry, at lowest order, is U(iV) x 
U(l)', with the additional condition that all gauge couplings are equal (to the SU(iV) 
gauge coupling). 

In the previous subsection, we introduced linear combinations of the extra gauge 
fields, A + and A-. Correspondingly, we may identify gauge groups U(l)+ and U(l)_ 
which are orthogonal combinations of the original extra U(l) and U(l)' gauge groups. 
These gauge invariances are described by the Ward identities 

^ + ^(4 + ^) = -^ (310) 

*lk = -7" B - (3 ' n) 

which follow from the Ward identity (3.3) and the phantom equation (3.8). 
The corresponding gauge-fixing conditions are 

^ = ™(d-A T + ZNB T ). (3.12) 
oB± g 2 

Note that the matter fields are neutral under U(l)+. This is in line with the obser- 
vation that A + does not interact at all. 

We will demonstrate in section 4 that interactions of A + are forbidden in the 
effective action to all higher orders, by the phantom equation, and there is a U(l)+ 
symmetry. On the other hand, regarding U(l)_, there is no guarantee that this 
symmetry can be enforced to all orders. If the gauge representation of the matter 
fields is chiral (QCD happens to be vector-like), the extra U(l)_ symmetry might 
be anomalous even if there is no anomaly of the original SU(iV). However, by the 
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argument in the previous subsection, we see that a U(l)_ anomaly will not contribute 
to observables, so it does not invalidate unitarity of the S matrix projected onto 
physical states. 

4 Renormalization 

In this section we outline a proof that the color-flow version of QCD can be extended 
consistently to all orders in perturbation theory. We follow the algebraic renormal- 
ization procedure, as explained in detail in Ref. [26]. Algebraic renormalization is 
mainly useful for proving the perturbative renormalizability to all orders of a class 
of models, in particularly for demonstrating the absence of anomalies, i. e. that sym- 
metries can be maintained in the quantized theory. 

It is convenient in practical loop calculations to employ a regularization that 
maintains the rigid and gauged symmetries of a Lagrangian in order to reduce the 
number of required counterterms. For recursive proofs to all orders, on the other 
hand, it is important to prove the absence of divergencies that would require the 
introduction of counterterms which would break a symmetry of the n-loop effective 
action. For this purpose one can use a non-invariant regulator for which the model- 
independent quantum action principles [27-30] have been established, which limit 
the possible counterterms to a finite set that can be enumerated explicitly. 

The inductive proof of renormalizability proceeds then in three stages: first 
one derives a set of Slavnov- Taylor (ST) identities for the effective action from the 
rigid and gauged symmetries of the model under consideration (cf. sect. 4.1). Subse- 
quently, one proves that the tree- level effective action (i.e., the Lagrangian integrated 
over space-time) is the unique solution of these functional equations (cf. sect. 4.2), 
up to field and coupling constant renormalizations. This establishes an induction 
hypothesis. For the induction step at n-loop order, one shows recursively that these 
ST identities can be maintained in order n + 1 of perturbation theory through a suit- 
able choice of local counterterms if they are satisfied by the n-loop effective action 
(cf. sect. 4.3). 

The induction step consists of several parts. One first restricts the terms that 
can possibly violate the ST identities at order n + 1, assuming that they hold at 
order n. The quantum action principles guarantee that these terms take the form 
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of local operators with well-defined dimension and quantum numbers, and the ST 
identities impose further algebraic constraints. The coefficients of the local operators 
can be computed as the values of certain Feynman graphs. If there are operators with 
nonvanishing coefficients, one has to check whether they can be cancelled by adding 
suitable non-invariant counterterms to the Lagrangian. This is a purely algebraic 
problem. If there is a solution, the ST identities can be restored by the counterterms 
and renormalizability is proved. Otherwise, the symmetry is manifestly broken by 
an anomaly. 

4.1 Conditions Imposed on the Effective Action 

First, we need a precise (perturbative) definition of the quantum field theory. We 
define the effective action of color-flow QCD as a solution of the usual ST identity 
for SU(iV), here extended to U(JV): 
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(d-A + £B). 



(4.2) 



with the linear gauge-fixing equation 

5T _ 2 
6B ~ f 

These identities involve the gauge, gauge-fixing, and ghost fields Ap B^, c*-, 3, ex- 
pressed in the color-flow basis, the matter fields and the sources for the BRST 
variations, p A , p 1 p, p^p, p c - 

Furthermore, the effective action depends on the phantom field A and its asso- 
ciated NL field B. The dependence on those fields is accounted for by the phantom 
equation (3.8), which we state here in the form 
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and by the gauge-fixing condition for A, 
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(4.3) 



(4.4) 



The algebraic conditions are supplemented by renormalization conditions which 
fix the field normalizations and a minimal set of interaction parameters. In the 
present model, these are just the gauge couplings. 
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Alternatively, in a minimal-subtraction scheme such as MS, we read off renormal- 
ization conditions after regularizing the loop integrals and subtracting the divergen- 
cies. This can also be done here, but with a caveat: we may use any renormalization 
scheme for the SU(iV) coupling and for the ordinary SU(iV) fields, but in order to 
keep the simplicity of the approach, we impose the condition that the singlet U(l) 
coupling is equal to the SU(iV) coupling and that the field normalizations are identi- 
cal. (The ST identity does not enforce this condition since U(iV) is not simple.) For 
the U(l)' coupling and the phantom field, everything is then fixed by the phantom 
equation and the phantom gauge-fixing condition. 

4.2 Lowest-Order Solution 

The next step in the renormalization procedure is to identify the lowest-order effec- 
tive action. We can verify that the Lagrangian (2.13) is the unique solution of the 
constraints of the preceding section (with correct gauge group and representation) 
which is an integrated local polynomial in the fields of dimension equal to four. 

The BRST invariance of the theory, and of the lowest-order approximation in 
particular, is encoded in the ST identity (4.1) and leads to lowest-order BRST trans- 
formations of the form 

sA = dc — i[A, c] sc = ice sc = B sB = 

sip = icip sip = iipc. (4-5) 

The U(l) part decouples. The corresponding ghost fields Trc and Trc are free and 
can be dropped, such that this part of the ST identity can be replaced by the U(l) 
Ward identity (3.3). 

4.3 Inductive Renormalization 

For an inductive proof of renormalizability, we can now assume that we have found a 
Lagrangian which generates a renormalized effective action at order n that satisfies 
all conditions. The renormalized order-n + 1 effective action does not necessarily 
have this property, yet. 

We do not repeat the standard derivation and just state that the SU(iV) ST 
identity (4.1) and gauge-fixing condition (4.2) can be extended to order n+1, possibly 
by adding non-invariant local counterterms to the Lagrangian. Clearly, we can handle 
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d"Tr —- + d v ^^ + — ( <9 2 Tr B - Nd 2 B) = d u Xl n+1) (4.7) 
5 A" r)A" a 2 V / 



the gauge- fixing condition (4.4) in analogy to the SU(iV) gauge-fixing condition, so 
this equation can also be established at order n + 1. 

The only non-standard part is the phantom equation (4.3). We do not know 
whether the right-hand side is zero at order n+1. The quantum action principles [27- 
30] tell us, however, that the equation must take the form 

ArO) Ar( n ) 9/ ~ \ 2 / ~\ 

Tr w + r TrG ^ - Nd " G ^y^ \ d ^ B - Nd » B ) = ^ n+1) > (4-6) 

where the effective action r^ n ^ is evaluated from the order- n Lagrangian, including 
terms in the result up to order n + 1, i. e., diagrams with n + 1 loops. The opera- 
tor xl n+1 ^ on the right-hand side is a local polynomial in fields and derivatives of 
dimension three with the same quantum numbers as the left-hand side. 
Taking the divergence of (4.6), we obtain 

This is precisely the Ward identity for the U(l)+ symmetry (3.11) with an anomaly 
on the right-hand side. We have to verify that the right-hand side is in fact zero. 

It is well known that the only possible obstruction to such a U(l) Ward identity 
is the Adler-Bell-Jackiw (ABJ) anomaly [13], i.e., an operator corresponding to a 
triangle diagram with three external gauge bosons. In a vector-like theory such as 
QCD, this vanishes. Nevertheless, we may consider a chiral SU(iV) theory where this 
is not necessarily the case. 

With the hypothesis that the phantom equation (4.3) is valid for the renormalized 
fields and interactions at order n, we can apply the change of basis (3.5) to the order-n 
renormalized fields and rewrite this equation as 

- ™ (fl^ - ff>Tr B W ) = 0, (4.8) 

in complete analogy with the tree-level equation (3.8). Taking the second derivative 
with respect to matter fields of this relation, we immediately see that A^ does not 
interact at all. In particular, the interaction (i/jA+'i/j)^ vanishes. Hence, there is no 
triangle diagram involving A + that can contribute to the ABJ anomaly. Integrating 
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this zero, we conclude that the right-hand side of (4.7) vanishes at order n + 1. The 
U(l) + symmetry is preserved. 8 

Since the right-hand side of (4.7) vanishes, we have d v Xv = 0. A divergence- 
less vector must be the derivative of an antisymmetric tensor, thus 

X (n+1) = dHYft+V. (4.9) 

Y( n+1 > is a SU(iV)-mvariant tensor of dimension two, a local polynomial of the fields. 
The only possibility is 

F> +1 ) = a^TrG^ + (4.10) 

with constants a^ n+1 ^ and b^ n+1 \ If we now modify the Lagrangian by the finite 
counterterms 

AC = -a^TrG^TrG^ + ~b {n+1) G^& u , (4.11) 

and re-evaluate Eq. (4.7), this term is cancelled. 

In summary, by adding local counterterms to the Lagrangian, we are able to 
satisfy all required conditions to order n+1. In particular, we establish the vanishing 
of all nontrivial interactions of A + , i.e., the exact cancellation of the singlet and 
phantom interactions. 

We note that the cancellation does not apply to the orthogonal combination A_, 
so in the presence of chiral matter, an anomaly operator involving is allowed. 
However, as argued in Sec. 3.2, as long as A + does not interact, an operator which 
involves A_, even if it formally breaks the U(l)_ gauge symmetry, has no observable 
effect and can be ignored. 

This completes the inductive proof of renormalizability. Unitarity of the S- 
matrix is then established to all orders by the ST identity and by the phantom 
equation, as argued in Sec. 3.2. 

In practice, adding those counterterms is necessary and natural. In the usual MS 
renormalization scheme, the singlet gluon propagator will receive a renormalization 
different from the octet gluon propagator, because only the octet has self-couplings. 

8 Actual calculations are carried out in the original basis (A, A). In this basis, the argument im- 
plies exact cancellation between the integrands of loop graphs involving A and A. This is analogous 
to the Standard Model, where the absence of anomalies is evident in the gauge basis (W , B), but 
involves a similar relation between 7 and Z interactions in the physical basis W^, Z,j. 



-19- 



In order to keep the color- flow scheme simple, we require finite counterterms that 
restore the equality of the propagator residues, and thus of the SU(iV) and U(l) 
couplings, at the next order. Then, we also have to renormalize the phantom propa- 
gator, i.e., the U(l)' coupling, to the very same value. This renormalization coincides 
with (4.11). We can do this freely since the values of the U(l) and U(l)' couplings 
do not enter any observable quantity. 

5 Applications: Tree-Level Amplitudes 

The color-flow approach is particularly useful for the automatic calculation of tree- 
level (squared) amplitudes. In this section, we choose a few simple examples that 
show how this works in practice. 

5.1 Algorithm 

For any amplitude, we replace the QCD diagrams with quarks and gluons by corre- 
sponding color-flow diagrams. These contain in place of each octet gluon, either a 
U(iV) gluon (double color line) or a phantom (no color). The latter appears only in 
places where the gluon directly connects two fermion (i.e., single-color) lines. Each 
vertex gets a factor l/\/2 due to the different normalization of the QCD coupling. 
After squaring the amplitude, the color lines of each contributing color-flow diagram 
are connected to the lines of the interfering complex-conjugated diagram, both for 
the incoming and the outgoing state. 

Instead of starting with ordinary QCD diagrams, an implementation may con- 
struct the amplitude diagrammatically from scratch, treating U(iV) gluons and quarks 
of definite color and phantoms as distinct, ordinary particle species. The cross sec- 
tion is computed by squaring diagrams, including all interference terms. We recall 
that U(iV) gluons do not interfere with phantom gluons. 

Each squared or interference diagram has a color weight W which is simply a 
signed integer power of N, given by 





where 



L = Number of distinct closed color lines 
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/ = Number of internal phantom propagators 
E = Number of external phantom particles 

The N L factor originates from contracting each string of Kronecker deltas, which 
represent color conservation at vertices, to a single 5\ = N. The —1/N factors in 
the phantom propagators, including those across the cuts, result directly from the 
Lagrangian. 

It is possible to absorb the (-1/N) 1 factor for internal lines by including — 1/N 
in the phantom propagator, as done in the Feynman rules in Appendix A. The other 
factors are applied after squaring the amplitude. Alternatively, we could absorb all 
1/N factors, but not their signs, in the phantom- fermion coupling. 

In automatic calculations, the color-flow approach has the advantage that the 
combinatorics of constructing Feynman graphs is already implemented, so the pro- 
liferation of diagrams does not raise a bookkeeping problem. (The algorithm should 
avoid to compute identical kinematics twice, however.) Counting distinct colors is 
rather simple and can be done, for tree graphs, by looking at the external state. In 
fact, common conventions [14] require an event-generating program to classify the 
external state in terms of color connections. On the other hand, computing color 
factors algebraically in a T a -f abc basis requires some additional infrastructure, and 
the transformation to the color-connection basis has to be done explicitly. 

Algorithms that do not expand an amplitude in diagrams but compute off-shell 
wave functions recursively [16-18] benefit even more from the color-flow approach. 
Assigning color factors to diagrams is not applicable there, so a straightforward 
implementation of QCD color algebra would require keeping a color degree of freedom 
in each off-shell wave function while constructing the amplitude, and applying a color- 
matrix multiplication at each vertex. By contrast, in the color- flow approach each 
color line can be understood as labeling an independent particle species, hence the 
algorithm need not know about color at all. It just has to distinguish particle species 
and their respective Feynman rules. 

In the following examples, we verify the color-flow result against the equivalent 
col or- algebra result, where we evaluate the appropriate trace of T a and f abc matrices 
directly. We are not interested in the kinematical part of the calculation that involves 
propagators, Lorentz factors, and integrations over momenta, so we omit them. We 
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also leave out the QCD coupling which is attached to each vertex, and the color- 
averaging factor for the initial state. We just quote color weights. 

5.2 Quark- Antiquark Scattering 

We are looking at the process (q ^ q') 

qq — > qq' (5.2) 

in pure QCD, so at tree level there is a single diagram, gluon exchange. 

To eliminate all open color indices, we square the amplitude and sum over colors. 
In the diagrams, squaring connects the final state of the amplitude and its complex 
conjugate; this is indicated by a vertical dashed line. Similarly, the initial states of 
the amplitude and its complex conjugate, the open lines at the left and right margins 
are also understood to be pairwise connected. 



i i 




i 



For this squared amplitude, standard color algebra yields 

W = Tr [T a T b ] Tr [T a T b ] = ^5 ab 5 ab = - (N 2 - l) (5.3) 
as its color weight. 

Color-flow QCD gives four diagrams for the squared amplitude, which differ only 
in the color factor. We recall that there is a factor N for each closed color line, and 
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a factor —1/N for each phantom propagator. Including l/y/2 for each vertex, we 
obtain 

w = j ( N2 + N (w) + N (tt) + ^ fir)') = i - L > ■ (5 ' 4) 

We observe that there is some redundancy in this case which could be eliminated 
before computing the result. In any case, we can read off the color factor directly 
without using color algebra. 

5.3 Four-Jet Production in e + e~ 

Let us now consider 

e + e" -> qqgg, (5.5) 

where we can ignore the colorless initial state when drawing color-flow diagrams. 
For the amplitude, we have three QCD diagram structures which decompose into 
ten color-flow diagrams: 




We now look at the squared color-flow diagrams. Symmetries between diagrams 
simplify the calculation. Summing over all diagrams with common kinematics, we 
recover the QCD color weights as linear combinations of powers of N. We show three 
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terms, the remaining ones are equivalent regarding their color flow: 




2x r=S=^ 
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In the color-algebra formalism, we compute the color weights 

rj-ij. ppaj"*] j arfc y fecrf = _ (j\r 2 i) 
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which agree with the color-flow results, as required. 



(5.9) 





(5.12) 
(5.13) 
(5.14) 



6 Applications: Effective Interactions From Loops 



At tree level, the color-flow method expresses all amplitudes in terms of Feynman 
rules that exclusively contain U(A r ) gluons and phantom gluons. This allows us to 
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compute the color factor of any color-flow diagram in the expansion of a squared tree- 
level amplitude by merely counting the number of distinct external color-flow lines, 
as detailed in the previous section. The factors for the internal phantom propagators 
are accounted for in the normalization of the propagators. 

Once loop diagrams are involved, this is no longer true. Obviously, there may be 
closed color loops not attached to external lines. As a further complication, there is a 
particular class of loop (sub-) diagrams where singlet gluons make their appearance, 
distinguished from the two classes of U(iV) and phantom gluons that we encounter at 
tree level. Fortunately, this only slightly complicates the algorithm, and computing 
color factors remains straightforward. 

We are particularly interested in loop amplitudes that can be inserted as effec- 
tive vertices in tree-level diagrams. We absorb the color sum over internal closed 
color lines in the corresponding vertex factor. As a result, we again can deduce the 
remaining overall color factor of any amplitude merely looking at the external lines. 

An important example is the coupling of Higgs and electroweak bosons to gluons, 
which occurs first at one-loop level. As in the previous section, we are interested only 
in color flow, so we ignore all kinematical and coupling factors. 

6.1 gg^H 

Let us consider colorless particles (in particular, the Higgs boson) coupled to a gluon 
pair via a fermion (more generically, matter) loop. From the color flows in the 
triangle diagram 




we can derive a set of equivalent "effective" Feynman rules 
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Thus, in the resulting one-loop effective action there appear vertices that explicitly 
couple phantom gluons to singlet gluons (and to themselves). In the effective action 
of pure QCD, this is not the case. Graphically, we identify the singlet, the projection 
of a U(iV) gluon, as a color line being reflected at the vertex. Note that there is no 
Feynman rule 



\ 



(6.3) 



because the singlet-singlet interaction is incorporated in the vertex where the color 
lines pass through. 

These new Feynman rules for the effective vertex yield the correct result for 
H — >■ gg and, equivalently, gg — > H. Squaring the amplitude, we get the set of 
squared diagrams (external Higgs lines not drawn): 
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Here, there are no interferences. To get the correct factors, we only have to remember 
that each phantom gluon comes with a factor —1/N, while the singlet gluon carries 
color, which yields a factor N when summed over. As in the tree-level case, we just 
have to count color lines crossing the cut to obtain the color weights of the squared 
diagrams. 

6.2 gg HH 

Things become interesting when there can be two loop-induced effective vertex in- 
sertions, as in the process gg —> HH. (The irreducible effective ggHH vertex has 
the same color structure as ggH, so we do not consider it here.) In this case, singlet 
gluons interfere with U(iV) gluons, projecting out the singlet part on the other side 
of the cut: 
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N 2 -l (6.5) 



iV^ J ' iV 2 ~ N 3 

Again, the correct factor is recovered, albeit many cancellations are involved. 

At this point, we digress somewhat and discuss a possible way to remove this 
redundancy before the sum is computed. Depending on the method by which the 
squared amplitude is computed, this may be useful for improving efficiency 

We would like to apply the direct cancellation between a phantom and a singlet 
line, where it appears obvious. If the amplitude is completely expanded in terms of 
diagrams, explicit singlets appear always attached to one of the effective vertices that 
involve a closed color loop, including the loop pertaining to the effective non-QCD 
interaction that we are considering here. 

Starting at this point (the vertex on the left in the figure), we could apply a 
procedure that looks like 



N 

(6.6) 

We can replace the color of an isolated singlet gluon that is attached to further gluon 
lines by a closed color loop, if we divide by 1/N. We then remove the closed color 
line from the diagram, attaching the factor N to the vertex on the left where the 
singlet originated. The same factor will be present in the analogous diagram where 
a phantom originates from the vertex. The two contributions now differ only in sign, 
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and thus cancel. 

This algorithm is nonlocal in nature and can be implemented, at face value, 
only if the amplitude is expanded in terms of diagrams. In automatic computational 
programs that construct the amplitude recursively, without expanding Feynman di- 
agrams, the cancellation procedure would have to be implemented as part of the 
recursive calculation. 

In any case, if we apply this argument to the gg — > HH amplitude before squar- 
ing, we obtain the much simpler result 




(6.7) 



When we compute the square, we have to take into account that singlets, being 
contained in the U(iV) gluon matrix, interfere with U(iV) gluons. This interference 
acts as a projection operator which turns a U(iV) gluon into a singlet, on the other 
side of the cut: 




Graphically, we can again apply the cancellation procedure and obtain the final 
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diagrammatic result 





AT. 



AT 



JV 



N 



N 2 - 1 (6.9) 



There is no redundancy left. Note that, effectively, the cancellation procedure has 
switched the sign of the phantom-loop graph, so we would have obtained this result 
if we just had included this graph with switched sign, but no graphs that explicitly 
involve singlets. This observation might be generalized and incorporated into the 
algorithm. 

6.3 H -)■ ggg 

The effective one-loop vertex H — > ggg has two color structures, f abc and d abc , which 
originate from the difference and sum of the possible loop orientations, respectively. 
In the color-flow basis, the Feynman rules are 
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The square of the / term is straightforward: 

TV 3 N 
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2iV(iV 2 - 1) (6.12) 
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The interference vanishes as expected. The square of the d term yields 
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(6.13) 



last two terms, i.e., 



Eliminating the external singlet lines as before would allow us to combine the 

27V - 127V + — = 2 ^iV - 4 

which is slightly simpler. 



27V 3 + 27V - 127V + — = 2 f TV - ) (TV 2 - 1) 



(6.14) 



7 Conclusions 

The expansion of QCD amplitudes in a color-flow basis has been known as a useful 
device in various contexts of perturbative and non-perturbative calculations. In this 
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paper, we have demonstrated that it can also be understood as a field theory of 
its own, a priori different from, but equivalent to, standard QCD. We have shown 
that this field theory is well defined, renormalizable, unitary, and yields observable 
predictions identical to ordinary QCD. 

For practical applications, color-flow QCD has advantages in the context of au- 
tomatic calculations, where it systematically generates color-connected amplitudes 
that can be matched to parton-shower and hadronization algorithms, avoiding some 
extra effort that pertains to color algebra and change of bases. In particular, it 
is useful for algorithms which do not work with a diagrammatic expansion of the 
amplitude. 

We have implemented color-flow QCD in the O'Mega matrix element generator, 
which is the tool for tree-level amplitude generation contained in the WHIZARD 
event generator package [19] together with its parton shower generator [32]. The 
color-flow approach also provides a convenient way of incorporating exotic color 
interactions. If desired, the methods presented in the present paper can readily be 
extended to other exotics, e.g., color-decuplet fields. 

For higher-order calculations, color-flow QCD can be taken at face value, pro- 
vided the renormalization procedure is properly implemented (cf. Sec. 4), which is 
essentially trivial. Concerning mixed QCD-electroweak processes at loop level, the 
examples in Sec. 6 show that there are minor technical complications, which however 
do not invalidate the procedure. A comprehensive treatment of NLO calculations 
would also require a color-flow exposition of (dipole) subtraction and parton split- 
ting kernels. This is not discussed in the present paper, but straightforward and part 
of a different publication [33]. 

Finally, we remark that the methods developed in the present paper could be 
generalized and applied to other gauge groups, including exceptional groups, using 
the completeness relations from [3]. 
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A Feynman Rules for Color-Flow QCD 



Here, we list the Feynman rules for SU(iV) gluons interacting with fermions in the 
fundamental representation, in the color-flow formalism. We shift the coupling g / y/2 
from the kinetic gluon terms in (2.2) and (2.13), respectively, to the vertices by 
renormalizing gluon and phantom fields accordingly, to obtain canonically normalized 
kinetic terms. Expanding in terms of components, the color-flow Lagrangian becomes 



+ B)(d ■ A\) - NB(d ■ A) + ^-B)Bl - N^-B 2 + C 



ghost 



(A.l) 



where 



(A.2) 



Note that in the color- flow formalism, the basic coupling emerges as g/y2. In the 
Feynman rules below, Lorentz and momentum factors are omitted, they retain their 
usual form. 

Propagators: 



a \£QQQQQQQQQQ£u b 

_ iS ab 



(A.3) 



_ iS ab 
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(A.4) 



Vertices: 
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kl mn kl mn 




op ji op 



kl mn kl mn 




op ji op 



kl mn kl mn 




ji op ji op 



B Color-Sextet Particles 

QCD and the Standard Model contain colored fields only in the fundamental and 
in the adjoint representation of SU(3). However, various Standard Model extensions 
that are currently under discussion provide extra colored fields and interactions. 
Automatic computation programs should be able to deal with such extensions, and 
the color-flow approach provides a straightforward means to achieve this. As a 
particular example, we describe the color- flow Feynman rules for color-sextet scalar 
particles. 
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Color-sextet particles or, more generally, particles in the symmetric bi-fundamen- 
tal representation of SU(iV) might exist at mass scales accessible at the LHC. A 
possible source are extended Higgs (super-) mult iplets in unified gauge theories of the 
strong and electroweak interactions. 

We consider a color-sextet scalar particle denoted as o. An explicit representa- 
tion is a complex symmetric matrix in color space with six independent entries 0^2,3 
and o-i 2,3, 



o 



( ^ 75^ 73^ 



(B.l) 



with components Oy, where o^ = Ojj. We also need the antiparticle o ij . This particle 
couples to gluons and gluon pairs. It may also couple to colorless states (e. g., Higgs), 
and linearly to quark pairs. The Lagrangian is 

C = D^DW - m 2 a t ^ - \Uo u cr il - g^q^f + h.c.) (B.2) 

The covariant derivative, in color-flow QCD, is given by 

D^o = d^a - i(^o + oAl - 2i M o) (B.3) 

We derive the Feynman rules. The propagator has two components which sym- 
metrize the color flow: 



i/ 2 xkxl i/2 xlxk 



(B.4) 



The same symmetrization is needed for external sextet particles. 

Gluons interact symmetrically with both color lines. The single-gluon interaction 
splits into three distinct Feynman rules: 



•ffh'rt 



1 



(B.5) 
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The double-gluon interaction splits into six distinct Feynman rules: 



SI 





W9tiv iff g^v 2i£ fg^y 



x6$8?6tf l p ^5 l 3 5°J k p *5°5yS l J k p 



V 




Couplings to colorless particles are simply 



i\S k S l 3 



(B.6) 



(B-7) 



and the coupling to quark pairs has the Feynman rules 

i 




19qq' S l S i 



(B.8) 
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